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Point to point correlators of currents are computed in quenched QCD using a chiral lattice fermion 
action, the overlap action. I compare correlators made of exact quark propagators with correlators 
restricted to low (less than 500 MeV) eigenvalue eigenmodes of the Dirac operator. In many cases 
they show qualitative resemblence (typically at small values of the quark mass and distances larger 
than 0.4 fm) and they differ qualitatively at larger quark masses or at very short distance. Lattice 
results are in qualitative agreement (and in the difference of vector and axial vector channels, 
quantitative agreement) with the expectations of instanton liquid models. The scalar channel shows 
the effects of a quenched finite volume zero mode artifact, a negative correlator. 



I. INTRODUCTION 



One of the earliest attempts to study nonperturbative behavior in QCD was through an analysis of correlators of 
^ [ hadronic currents. The QCD sum rule approach treats short distance physics in the context of the operator product 
^ ■ expansion (OPE) and parameterizes long distance physics in terms of vacuum condensates 0,1). The approach was 
I ' quite successful in the vector and axial vector channels, but nonperturbative effects are very large in the pseudoscalar 
. and scalar channels. The realization that instanton effects appear in precisely those channels in which nonperturbative 
' effects are large led to the development of instanton liquid models ||^,^| . It would be an interesting exercise to "validate" 
, the parameterization and results of these models directly from the QCD Lagrangian, by using lattice methods. This 
' paper is a step in that direction. 

We are concerned with flavor non-singlet current correlators separated by a (four dimensional Euclidean) distance 

X, 



Q-i' 

^ where the current will be proportional to 



U,{x) = Tr(Jf (x)Jf(O)) (1) 



J^x) = i,{x)T^mi^{x) (2) 



K> , with T° an isospin label and T a product of Dirac matrices. In practice (and to be consistent with the standard 
' approach to these correlation functions Q) I present results for the correlators scaled with respect to the massless 
[ free field current-current correlator 11° (x) (the precise kind of free fermions will be described below) 

R,{x)^U,{x)/U°{x). (3) 

The label i will include the pseudoscalar (PS), scalar (S), vector (V) and axial vector (A) currents. 

The Euclidean space correlators whose Fourier transforms into coordinate space I am measuring obey a dispersion 
relation which relates them to their absorptive part in Minkowski space, the spectral density p{s), 

U(,^)^'- [dsJ^^. (4) 
TT J s + q-' 

p{s) is in turn is proportional to the total cross section for scattering in the appropriate channel, and should then be 
positive. The point-to-point correlators measured here are then related to the spectral density through 



n(x) = / d^qexpUq ■ x)— [ ds \ (5) 
TT J s + q^ 



or, more directly. 
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n(x) = i y dsp{s)D{^s, \x\) (6) 

where 

^(v^,N) = A r^=exp(-Nv/^) (7) 

is the EucUdean propagator for a free particle of squared mass s. If we were measuring a conventional (for the lattice) 
g = correlator, Eq. ^ would be replaced by 

2m 

In principle, lattice calculations could check results of instanton liquid models Q. In practice, however, lattice 
simulations have been contaminated by the use of fermion actions which are not chiral. Since the continuum analyses 
depend crucially on chiral symmetry arguments, one might wish to be cautious about drawing conclusions from 
simulations with nonchiral lattice fermions. The recent discovery of lattice actions which support an exact chiral 
symmetry Q (notably the overlap action 0) allows one to revisit these questions in a theoretically clean context. 
That is the subject of this paper: a study of current correlators from an overlap fermion action. 

Because I am working with an action with exact chiral symmetry, I have to be careful with my definitions and 
make them consistent with current commutation relations. If I assume that I am modeling a theory with a flavor 
SU{2)v X SU{2)a — 0(4) symmetry, I can write down two real 0(4) vectors of currents 

'?l = (^^/o) (9) 

and 

$2^{a^o,v) (10) 

where 7r° = ii!^]^{T°-)4' = Jps is the pseudoscalar (PS) current, ag — ~{1>{t°- /2)%Ij = Js is the scalar (S) isotriplet 
current, /o = V'V' is the scalar isoscalar current, and 77 = iipj^ip is the pseudoscalar isoscalar current. (Correlators of 
the last two currents involve disconnected diagrams and I will not consider them further.) 

In the cases of the V (vector) (F — -fi^i), and A (axial vector) (F = jfj,^5) currents, I sum the correlator over the 
four values of fj,. The most interesting way to present results for correlators in the vector and axial vector channels is 
to look at the sum or difference of vector and axial vector correlators 

Uv{x)±UAix) 

'^^^^("^ ^ 2n°.(.) ■ (^'^ 

In the sum rule/OPE approach Ry+A is dominated by perturbative physics and is expected to take a value very close 
to unity, while Ry-A is zero at small x and receives only nonperturbative contributions which are relevant to chiral 
symmetry breaking. The instanton liquid model produces a large part of the correlator Rv~a, as emphasized recently 
by Schafer and Shuryak [pL 

In a recent publication 1 3] we have shown that the low lying eigenmodes of a chiral lattice fermion Dirac operator, 
an overlap fermion operator (described in Ref. |ic|] ), has a local chiral density 4!{xyj5^p{x) which shows a peaked 
structure. The positions and signs of the peaks are strongly correlated with the locations of topological objects, which 
would be identified as instantons and anti-instantons detected using a pure gauge operator. Zero modes correlate 
with only one sign of topological objects, while nonzero eigenmodes of the Dirac operator interpolate between both 
signs of topological object. This correlation dies away slowly as the eigenvalue of the mode rises. Spatially averaged 
correlation functions of hadrons made of light quarks are saturated by propagators of quarks restricted to a few low 
eigenmodes. The whole picture is very reminiscent of an instanton liquid model. 

Ref. did not consider the possibility that the density of the fermionic modes might be large in places where the 
local chiral density was small. Subsequently, this point was raised by the authors of Ref. . Ref. showed that 
this possibility did in fact not occur for the overlap action used in Ref. Q . Similar results have been presented using 
an alternative chiral lattice action (the Wilson overlap action) jl^ two lattice actions with improved but inexact 
chiral symmetry (domain wall fermions and an approximate Ginsparg- Wilson action |l5[| ) and with an improved 
operator in a lattice action with inexact chiral symmetry (the clover action) ||l6|. None of these last three works 
compared the distribution of chirality with topological charge density measured with a gauge observable. 
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It is an uncontrolled approximation to replace a quark propagator by the quark propagator restricted to a sum 
over a small number of small eigenvalue eigenmodes of the Dirac operator. However, if an observable computed with 
truncated propagators resembles the same observable computed with exact propagators, and if we further have seen 
that the low modes couple to particular structures in the QCD vacuum, then this resemblence is a strong (qualitative) 
signal that the particular vacuum structure is connected with that observable. 

For example, consider the difference of vector and axial vector currents. Rv~a{x), either extracted from r— decay or 
measured in Monte Carlo simulation, is zero at small x and increases with x. The same behavior is seen in a calculation 
based on propagators truncated to a small number of eigenmodes: in fact, at small quark mass (pseudoscalar to vector 
mass ratio less than 0.5) the difference between the exact and truncated Ry-Aix) is small. Since in |9| the locations 
of peaks in the low eigenvalue modes were correlated with the locations of peaks in the topological charge density, 
one can infer that instantons are connected with the rise in i?v'_^(a;). 

If the current correlator based on truncated propagators does not reproduce the full calculation, the truncated 
propagator must be missing some important physics. For example, in most channels, R{x) approaches unity as x falls 
to zero. Approximating quark propagators by a few low eigenvalue modes does not reproduce this behavior: in that 
case R{x) typically vanishes at small x. Presumably R{x) — > 1 as a; — *■ is asymptotic freedom at work: the correlator 
is reducing to its free-field value. The low eigenmodes are extended in space and decouple from short distance physics. 

Previous studies [ pT|Jl8| of point-to-point correlators have mostly been concerned with the long-distance behavior 
of these observables, including the extraction of particle masses and couplings. I will not attempt to do that here, 
because my simulation volume is very small. Some studies have also compared the correlators in a "full" QCD 
simulation with propagators truncated to a set of low-eigenvalue eigenmodes. They see results similar to the ones I 
report, the main difference being that they need many more eigenmodes to saturate the correlator. I believe that this 
difference is due to the poor chiral behavior of the action, the usual thin link Wilson action, used by these authors. 

Along the way, I observe an interesting artifact of the quenched approximation in the scalar channel-the correlation 
function becomes negative. This behavior is incompatible with a normal spectral representation. In my data it seems 
to be associated with a finite- volume quenched approximation artifact arising from the exact zero modes of the Dirac 
operator. 



II. THE LATTICE CALCULATION 

The overlap action used in these studies is built from an action with nearest and next-nearest neighbor couplings, 
and APE-blocked links |^0|. Eigenmodes of the massless overlap Dirac operator D{0) are constructed from eigenmodes 
of the Hermitian Dirac operator H{0) = 75_D(0), using an adaptation of a Conjugate Gradient algorithm of Bunk 
et. al. and Kalkreuter and Simma |2^]. These eigenmodes are used to precondition the calculation of the quark 
propagator (by use of a conjugate gradient algorithm) and are additionally used to construct quark propagators 
truncated to some number of low lying eigenmodes. 

Later, we will need a few simple facts about the eigenmodes of the overlap Dirac operator, so let us recall them 
now. The eigenmodes of any massless overlap operator are located on a circle in the complex plane of radius xq 
with a center at the point (a;o,0). The corresponding eigenfunctions are either chiral (for the eigenmodes with real 
eigenvalues located at A = or A = 2xq) or nonchiral and paired; the two eigenvalues of the nonchiral modes are 
complex conjugates. The massive overlap Dirac operator is conventionally defined to be 

D{m) = {I )D{0) + m (12) 

2x0 

and it is also conventional to define the propagator so that the chiral modes at A = 2xq are projected out, 

D-\m) = l^-^(D-\m) - (13) 

^ ' l-m/{2xo) ^ ' 2xq' ^ ' 

Then the contribution to the propagator of a single (positive chirality) zero mode in the basis where 75 = diag(l, — 1) 
is 

I)- (14) 



The jth pair of nonchiral modes contributes a term 
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where, defining /i = m/{2xo), ej = Xj/{2xo), the entries are 



1 m(i - s') 

2X0 e] + y?{l^e]) 



(16) 



^^"-2xoe? + ;,2(i_,2) (17) 
and D{0)'^(pj — X'^cfij] the eigenmodes of -0(0) have eigenvalues 2xo{e'^ ± ~ e^)- For a summary of these useful 



(--J — hii<o ci5ciiiiiwvxt;o wi J-^ \^ ) iia-vc eigenvalues 2xo(cj -1- f^^jyj -L c j y 

formulas, see Ref. |Q (for the special case Xq — 1/2). 

The data set used in this analysis are generated in quenched approximation using the Wilson gauge action at a 
coupling P = 5.9. The nominal lattice spacing is a = 0.11 (inferred from the Sommer parameter using the interpolation 
formula of Ref. [^) or 0.13 fm from the rho mass. This will be discussed in more detail below. It consists of 20 12"*^ 
configurations. The fermions have periodic boundary conditions in the spatial directions and anti-periodic temporal 
boundary conditions. I calculated the ten smallest eigenvalue modes of H^{0) in the chiral sector of the minimum 
eigenvalue, and reconstructed the degenerate eigenstate of opposite chirality of H^{0) when one was present. These 
modes are then recoupled into eigenmodes of -D(O). Their eigenvalues have imaginary parts ranging up to 0.3/a-0.35/a, 
or about 500 MeV g. 

The lattice analysis has one strong point and one weak point. The strong point is the use of an action with exact 
chiral symmetry. This means that there is no additive renormalization of the quark mass. It also means that there 
are no exceptional configurations, so the simulations can be performed at small quark mass. Instead of exceptional 
configurations, there are contributions to the quark propagator from zero modes of the Dirac operator. Because of 
chiral symmetry, these zero modes contribute in a "dangerous" way, as finite volume effects only in the PS and S 
channels. The particular choice of action, with a fat link, also has rather mild multiplicative renormalization of current 
matrix elements. The perturbation theory has not been done for this action, but similar calculations done for the fat 
link clover action have vector and axial vector lattice-to-continuum renormalization factors quite close to unity . 

The weak point of the calculation is its use of a small volume. The lattice has a size of about 1.5 fm, assuming a 
lattice spacing of 0.13 fm. 

One annoyance encountered in this study is the different finite volume effects suffered by the free massless current 
correlator and the correlator in the nontrivial gauge background. At large distance the free correlator can receive 
contributions in which the quark and antiquark wind in opposite directions around the lattice, so that the correlator 
has the topology of a single line encircling the simulation volume. This contribution is not present in the nontrivial 
background, since it corresponds to a single Polyakov loop winding around the lattice. Confinement forces the 
expectation value of this operator to be zero. The way I got around this problem is the same as was used in Ref. [ pT| : 

Think of the finite volume lattice as a piece of an infinite lattice. A propagator from a point source in the finite 
lattice corresponds, on the infinite lattice, to a sum of propagators from sources which include the original source, 
plus a set of image points on all the "copies" of the finite lattice which tile the infinite lattice. We only want to include 
contributions to correlators from infinite- volume propagators where the source points for the quark propagator and 
antiquark propagator coincide, and so do the sink points. These quark-antiquark correlators look like closed fermion 
loops. Contributions in the finite volume in which the quark and antiquark propagate in opposite directions around 
the finite volume to get to the same sink point are topologically equivalent to lines in the infinite volume, where the 
quark source and antiquark source sit on different image points and the two propagators terminate on the same sink 
point. 

To compute the necessary free field current-current correlator, I construct an approximation to the infinite volume 
free fermion propagator by computing the free (overlap) fermion propagator on a large lattice (I used a 24"* volume) 
and approximate the quark-antiquark correlators in the small volume as a sum of the "direct" term and propagators 
from sources located on (the same) nearest image point. 

The hadron correlators also suffer finite volume effects. In particular, the lattice volume cuts off the long distance 
part of the hadron correlator, which in most cases contains the contribution to the correlator of the lightest state in 
the channel. Thus, I have not tried to do any fits to the mass of this state, in contrast to Refs. ||l^,|l^. Note however, 
that for short distances (less than 0.5 fm or so), the effects of finite volume are not too important. Also, comparisons 
of a full correlator to a correlator computed using truncated propagators are not affected by the free field correlators, 
so the observation that low modes saturate a correlator (or not) are not so compromised. 

At very light quark mass, chiral symmetry breaking begins to be modified by the small volume. The relevant 
parameter is the Leutwyler-Smilga Eq] parameter, z = nig'SV for quark mass to^, infinite volume condensate S = (qq), 
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and simulation volume V. A large volume corresponds to 2; 1. Using the value of S computed for this action at 
this coupling in Ref. [|l0|, z ~ 0.8(am^/0.01)) (to show z in units of the smallest quark mass used). At this quark 
mass the scaled lattice size L is itit^L ~ 2.4. 



III. RESULTS 



A. The pseudoscalar and scalar channels 

In Fig. |l| is shown the pseudoscalar correlator for several choices of light quark masses, with pseudoscalar to 
vector meson mass ratios (mps/mv) inferred from Ref. (a) aruq — 0.01 {mps/mv — 0.34); (b) aruq = 0.02 
i'mps/'n^v — 0.50); (c) airiq — 0.04 [mps/mv — 0.61); (d) arUq — 0.06 [mps/mv — 0.64). The separation x is 
measured in units of the lattice spacing in the simulation. Also shown is the contribution to the correlator from the 
lowest ten modes of the Dirac operator, as well as the contribution to the correlator given by the zero modes alone. 
Quark interactions in this channel are strongly attractive. The enormous attraction in the pseudoscalar channel is well 



known in continuum (compare Ref. [p6| ) and lattice (compare Ref. |17|) analyses. In instanton liquid phenomenology 
it is ascribed to the effect of the 't Hooft interaction, which is strong and attractive in the pseudoscalar channel. At 
low quark masses {mpv/mv < 0.5 or so) the low mode truncation saturates the correlator at larger x. As these are 
the modes which couple to instantons, there seems to be a connection between instantons and the strong attraction 
in the pseudoscalar channel at x/a = 5 — 9. 

The fact that the contribution of the low modes begins at R{x = 0) = and not 1 is just a sign that at tiny x 
the channel is dominated by free- field modes (which presumably do not couple to topology), but by a distance of 2-4 
lattice spacings (0.2-0.4 fm) the light quark mass correlator is completely dominated by the low modes. As the quark 
mass falls to zero, the pure zero modes make an ever larger contribution to the point-to-point correlators, as well: 
their contribution scales a 1/m^ times a quark- mass- independent function of x and y, 

(:r,y) = ^| ^ 0,(x)^,(y)tp. (18) 
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FIG. 1. Saturation of the point-to-point pseudoscalar correlator by low-lying eigenmodes of //(O)^. (a) am, — 0.01 
(mps/mv ~ 0.34); (b) aniq = 0.02 {mps/mv ~ 0.50); (c) am, = 0.04 {mps/mv ^ 0.61); (d) am, = 0.06 (mps/mv ~ 0.64). 
Octagons show the full hadron correlator. Squares show the contribution from the lowest 10 modes. Diamonds show the 
contribution from the zero modes, which just scales as l/(amq)^). 

As the quark mass rises, the approximation of Rps{x) by a small number (ten) of low- lying eigenmodes becomes 
poorer and poorer. This is shown in Fig. ^. The correlators flatten and approach unity over a wide range of x. That 
the low eigenvalue modes do not saturate the correlator shows physics of these heavier pseudoscalars has little to do 
with fermion eigenmodes which are strongly coupled to instantons. This is not a surprise. It points up the dangerous 
possibility, however, that lattice simulations involving pseudoscalars performed at a pseudoscalar to vector ratio above 
0.7 or so might give misleading results when extrapolated to small quark mass. 
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FIG. 2. Non-saturation of the point-to-point pseudoscalar correlator for heavier quark masses by low-lying eigenmodes 
of H(Of. (a) am, = 0.10 (mps/mv ~ 0.75); (b) am, = 0.15 (mps/m^ 0.84); (c) am, = 0.20 (mps/mv ~ 0.87); (d) 
airiq = 0.25 {mps/mv — 0.91). Octagons show the full hadron correlator. Squares show the contribution from the lowest 10 
modes. 



The scalar correlator shows similar qualitative behavior to the pseudoscalar correlator: strong deviation from unity 
away from small x and at small quark mass a large contribution from zero modes. See Fig. ^j. 
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FIG. 3. Saturation of the point-to-point scalar correlator by low-lying eigenmodes of H{0)^. (a) arUq = 0.01 
(mps/mv ^ 0.34); (b) am, = 0.02 [mps/mv ^ 0.50); (c) am, = 0.04 {mps/mv ^ 0.61); (d) arUq = 0.06 {mps/mv ^ 0.64). 
Octagons show the full hadron correlator. Squares show the contribution from the lowest 10 modes. Diamonds show the 
contribution from the zero modes, which again just scales as l/{amq)^ . 



Instanton liquid models predict a repulsive interaction in this channel. We see that, indeed, the interaction between 
quarks is strongly repulsive - so strongly repulsive, that at small quark mass, the correlator becomes negative at 
larger x. 

The bulk of this effect is due to the zero modes. To make the argument, I begin by considering the susceptibilities, or 
spatial integrals of the correlators [p2[ : In the pseudoscalar channel, configuration by configuration, if the configuration 
has topological charge Q, 



x,y,a 



V 



75 



(19) 
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(20) 



(21) 



The last line is the Gell-Mann, Oakes, Renner relation, which is satisfied exactly by overlap fermions. The isotriplet 
scalar susceptibility is 



8 



Xao - E (Jsi^)Jsiy)) = -|Tr (22) 

= -^^ + V^(^^--^)- (23) 

Note the negative contribution of the zero modes to the scalar channel. In general, one would expect the sum over the 
nonzero modes to contribute a result of order unity to the susceptibility. The contribution of the zero modes should 
scale away as 1/VV, since the existence of a topological susceptibility means that (Q^) ~ V. Thus the zero mode 
contribution to both the pseudoscalar and scalar channels is a finite volume artifact, but one which dies away quite 
slowly. 

So much for the integral of the correlator; what about its value point by point? The contributions where both of 
the quark lines propagate through zero modes are equal in magnitude and opposite in sign in the pseudoscalar and 
scalar channels. If the correlator is dominated by modes whose eigenvalues are large compared to the quark mass, one 
would expect the nonzero mode sum to contribute positively to the correlator. This is certainly the case for the free 
theory, and so we naturally have a situation where the free scalar correlator is positive and the interacting correlator 
is dominated by its zero modes and is negative. 

One can check that the zero mode is the source of the loss of spectral positivity by considering the scalar propagator 
in the subset of our lattices which have zero topological charge (5 of the 20 lattices). This is shown in Fig. ^. Notice 
that the scalar correlator in this sector is still saturated by low eigenmodes of the Dirac operator for mpg/my < 0.5 
at a distance of 5 lattice spacings {x ~ 0.55 fm). In this sector the channel only becomes strongly repulsive at heavier 
quark masses, where the low modes do not contribute to the correlator. 

The Q = sector of QCD is not QCD and so it remains an open question, what will happen at larger volumes. 
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FIG. 4. Comparison of the point-to-point scalar correlator to the scalar correlator composed of low-eigenvalue quark modes 
in the (5 = sector, (a) am, = 0.01 {mps/mv — 0.34); (b) aniq = 0.02 (mps/mv — 0.50); (c) aniq = 0.04 {mps/mv — 0.61); 
(d) aniq = 0.06 {mps/mv — 0.64). Octagons show the full hadron correlator. Squares show the contribution from the lowest 
10 modes. 



Finally, one might worry that all of the attraction in the pseudoscalar channel is due to the zero modes. The sum 
of pseudoscalar and scalar correlators, 



Rps+s{x) 



(24) 



has no zero mode contribution. 

This quantity is shown in Fig. ||. The low modes do not contribute to the correlator at small x. At larger x they 
show a significant attractive interaction, saturating the full-propagator correlator at larger x for mps/mv < 0.5. At 
heavier masses they produce too much attraction at larger x. 
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FIG. 5. Saturation of the sum of point-to-point pseudoscalar and scalar correlators sector by low-lying eigenmodes of H{0)^. 
(a) arriq = 0.01 {mps/mv ^ 0.34); (b) am, = 0.02 (mps/mv 0.50); (c) amq — 0.04 {mps/mv ^ 0.61); (d) amq — 0.06 
{mps/mv — 0.64). Octagons show the full hadron correlator. Squares show the contribution from the lowest 10 modes. The 
curve in panel (a) is from a simple phenomenological parameterization of the spectral density. 

Other possibilities for a negative spectral weight could exist in the scalar channel. Schafer and Shuryak |^ actually 
show that the scalar correlator in interacting instanton models becomes negative at a; ~ 0.5 fm. They ascribe this 
behavior to the strong repulsion in the scalar channel expected in instanton models becoming a little too strong, and 
their calculations of interacting instantons which include dynamical fermion effects soften the interaction and restore 
positivity. 

Recently, a discussion of the same effect has been given by Thacker et al |^ in the context of quenched chiral 
perturbation theory In full QCD the quark and antiquarks in the eta-prime, the flavor singlet pseudoscalar 

meson, can annihilate and the (multi) qq pair(s) Fock state of the meson can mix with a quarkless intermediate state. 
This mixing shifts the mass of the eta-prime. In the quenched approximation the sum of mixing graphs truncates 
with a single term and the propagator consists of a single pole and a double pole with a negative coefBcient, a ghost. 



G{p) 



1 



1 



1 



3 p^ +ni^ 



(25) 



A negative spectral weight in the scalar channel can appear because the scalar particle can couple to an intermediate 
state which is a combination of a pion and a quenched eta-prime, whose propagator is the second term (hairpin) part 
of Eq. E5l This is a connected graph so in the quenched approximation it can be constructed. It is shown in Fig. o. 
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(a) (b) 

FIG. 6. The quark line graph (a) and associated Feynman graph (b) which makes a ghost contribution to the scalar 
correlator. 

While we do not know the overall magnitude of the graph, the particle in all the propagators is the flavor nonsinglet 
pseudoscalar, whose mass is known from conventional spectroscopy calculations. If I assume that the scalar particle- 
pion-eta coupling isjust a scalar constant g, then I can compute p{s) from the graph of Fig.^(b). Because of the 
minus sign in Eq. ^ the amplitude is negative and it is also divergent at threshold, 

p{s) = ^ (26) 

^s(.s - 4m2) 

Then the shape of R{x) can be computed exactly up to an overall normalization. This contribution can have nothing 
to do with the zero modes because the spatial correlations of the zero modes do not depend on the quark mass. 
Attempts to fit the lattice data of Fig. ^ to a conventional [|| combination of continuum-like background p{s) = 
(3s/(87r2))0(s — So) plus a possible scalar resonance plus the ghost term (even varying the strength of the ghost 
coupling with quark mass) were unsuccessful. 

Schafer and Shuryak have presented predictions for point to point correlators in the instanton liquid model. It would 
be an interesting exercise to compare lattice results to theirs. Because of the presence of the zero mode contribution in 
the pseudoscalar and scalar channels, it probably makes the most sense to consider the sum of pseudoscalar and scalar 
channels, from which the zero modes decouple. Next, one must choose a lattice spacing. Because the correlators rise up 
so strongly with x, they are very sensitive to the choice of lattice spacing. I have collected a number of observables for 
this simulation, which happen to have more or less well-known continuum values, and which might be used to infer a 
lattice spacing, in Table |. The quantity tq is the Sommer parameter, from Ref. psf . The rho mass, pseudoscalar decay 
constant, and pseudoscalar matrix element A^r = {0\4'^5ijj\PS) (all extrapolated to zero quark mass) are determined 
from simulations on 12'^ x 24 lattices at the same parameter values while the infinite volume E measurement is 
from Ref. The phenomenological estimate for A^, which is equal to f-Kni^/ {niu + md), is from Ref. The 

inferred lattice spacing, shown in the fourth column of the Table, includes the unknown (but believed to be close to 
unity) lattice-to-continuum renormalization factors appropriate to each observable. 

The dominant feature of Rps{x) or Rps+s{x) is the pion. Its contribution is proportional to A^, or to the fourth 
power of the lattice spacing. This is the source of the sensitivity of the correlator to a. 



quantity 


lattice value 


continuum 


lattice spacing (fm) 


ro 


see II 


0.5 fm 


0.11 


mp 


amp = 0.50(2) 


770 MeV 


0.13 




af^ = 0.078(1) 


131 MeV 


Q.12/Za 


Att 


a'^K = 0.143(3) 


(470 MeV)^ 


0.16 ZaZttl 


(^^) 


a^E = 0.0039(1) 


(250 MeV)^ 


0.13 Zli"^ 



TABLE I. Table of observables and inferred lattice spacings. 
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I will somewhat arbitrarily take the lattice spacing to be 0.13 fm in Fig. |7|, where I compare the lattice results to 
instanton liquid model data from Ref. The solid line is a very naive phenomenological calculation of Rps+s'- I 
take a spectral function which has a delta-function pion of weight (with its value from the Table) and add in a 

step function representing a continuum: Pcontix) ~ 9{s — sq) with sq — 0.5 GeV^. This parameterization does a very 
good job of modeling the low quark mass lattice data, when one inputs the observed lattice values of pseudoscalar 
mass and A^, as is shown in panel (a) of Fig. In Fig. the data overshoots the model by about thirty per cent. It 
will take simulations at a smaller lattice spacing to disentangle scale violations from lattice spacing uncertainties, to 
improve this picture. 




0.0 0.2 0.4 0.6 0.8 

X, fm 

FIG. 7. Comparison of the sum of lattice point-to-point pseudoscalar and scalar correlators, extrapolated to zero quark 
mass (octagons), predictions of the instanton liquid model (diamonds) and a simple phenomenological model for p(s) (solid 
line). 



B. The vector and axial vector channels 



The potentially interesting physics in these channels involves a comparison of lattice data to the predictions of 
the OPE. With my sign conventions, the coordinate-space quark propagator (with flavor label a, b in slowly varying 
external fields is (for references see [^) 
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where the rather schematic "G/a:^" represents a long expression which will give the gluon condensate term in a 
mesonic correlator. With this quark propagator all the quark mass and quark condensate dependence in the "sum" 
Rv+A correlator cancels out, while the "difference" correlator for massive quarks is 



R 



V-A 



- imT'^{qq)x + Ax^ 



(28) 



and the Ax^ term is Ax^ — {tt^ /9)as{x) (qq)'^ x^ which is estimated by Schafer and Shuryak to be about (x/0.66 fm)^. 
The "sum" correlator sees the gluon condensate 



Rv+A — 1 



as{x) 



384 



27r= 
81 



as{x){qq) log(x )x'' + 



(29) 



Figs. H and || show Rv±a correlators from the full simulation and from the lowest ten eigenmodes. The low modes 
clearly saturate the "difference" correlator at mps/mv < 0.5. The "sum" correlator is essentially unity, and receives 
little contribution from the low eigenmodes. 



C 1 



C i 
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x/ a 
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FIG. 8. Saturation of the difference of point-to-point vector minus axial correlators by low-lying eigenmodes of H{0)^. 
(a) anig = 0.01 (mps/mv ^ 0.34); (b) am, = 0.02 (mps/mv ^ 0.50); (c) am, = 0.04 (mps/mv 0.61); (d) am, = 0.06 
(mps/mv — 0.64). Octagons show the full hadron correlator. Squares show the contribution from the lowest 10 modes. 
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2 4 6 8 10 C , 1 6 8 10 

x/a x/a 



FIG. 9. Comparison of the summed point-to-point vector and axial vector correlators from the full quark propagator and with 
the propagator built of low-lying eigenmodes of H{0)^. (a) am, = 0.01 {mps/mv — 0.34); (b) arUq — 0.02 {mps/mv — 0.50); 
(c) arriq — 0.04 [mps/mv ^ 0.61); (d) arrig = 0.06 {mps/mv ^ 0.64). Octagons show the full hadron correlator. Squares 
show the contribution from the lowest 10 modes. 
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show the lattice correlators with 
3y Ref. M (lines). In these figures 



Now the lattice data is collected at several values of quark mass, all of which are greater than the physical quark 
mass, and so to make connection with phenomenology it is necessary to extrapolate all the data to zero quark mass. 
The extrapolation itself can shed light on the OPE prediction: 1 will assume that point-by-point in x, R(x) = 
To -l- riijiq -f r2m^. For both correlators, the mass dependence is weak and the extrapolation is readily performed. 
However, the coefficients of vriq and expected from the OPE, Eq. are only seen at very low a;, x/a < A 

The extrapolated data can be compared to non-lattice results. Figs. |l^ and 
the instanton model of Ref. ||^ (crosses) and ALEPH r— lepton decay, as extracted 
I have re-introduced a physical distance scale (in fm) for the separation x by using the same lattice spacing (0.13 fm, 
from the rho mass) as used in the last section. Because the correlator is such a flat function of x, it is not too sensitive 
to the choice of lattice spacing. 

Lattice Monte Carlo data, the instanton liquid model, and the extracted tau decay data for Rv~a all agree nicely. 
This is a channel which has no short distance contributions, and the long distance contributions are dominated by 
instanton-sensitive eigenmodes. 

The lattice data for Rv+a also agree with the instanton liquid model out to a distance of about 0.6 fm, and then 
become less attractive. They undershoot the tau decay data by about five cent. It is plausible to assume that the 
latter discrepancy is just due to the lack of lattice physics below the cutoff scale, and to the fat link used in the gluon 
vertex, which smears the action out to about two lattice spacings (0.26 fm, nominally). It is the same discrepancy as 
the instanton liquid model data, which does not "naturally" include the perturbative Ofg/Tr part of the correlator. Of 
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course, as I have remarked earlier, I have not computed the lattice to continuum conversion factor for the vector and 
axial currents. While it is unlikely to account for all the discrepancy, it might be the source of some of it. 

0.3 I I I I \ \ I I I \ \ \ I I \ \ \ I I I I 




0.0 0.2 0.4 0.6 0.8 

X, fm 

FIG. 10. Comparison of the difference point-to-point vector and axial vector correlators from the overlap action (octagons), 
extrapolated to zero quark mass, and from the instanton model of Ref. (crosses) and ALEPH r— lepton decay, as extracted 
by Ref. [| (lines). 
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0.0 0.2 0.4 0.6 0.8 

X, fm 

FIG. 11. Comparison of the sum of point-to-point vector and axial vector correlators from the overlap action (octagons), 
extrapolated to zero quark mass, and from the instanton model of Ref. (crosses) and ALEPH r— lepton decay, as extracted 
by Ref. [| (lines). 

The authors of Ref. [|| perform fits to Rv±a{x) in order to extract the various condensate terms (which will appear 
as different powers of x"), as well as to look for logx^ terms so characteristic of the operator product expansion. I 
tried to do this with the lattice data. Because the lattice data and instanton liquid model data are so similar, the 
results will be similar to those found in Ref. ||^-that is, difficult to reconcile with the OPE. 

For example, consider Ry-A- The action itself has order discretization errors, which generate terms in R{x) 
proportional to x^. Good fits to the lattice data can be found beginning with any even power of x. However, if we 
take the continuum theory seriously and consider only fits beginning with x^ term, one can find fits for x < 0.4 fm or 
so of the form (x/ 0.6 fm)^, quite close to the expectation of the OPE and the results of Schafer and Shuryak | p7| |. 

Fits to Rv+A designed to expose the gluon condensate term are also not very successful. The fit parameters are 
quite sensitive to tie range of x. Most fits (which range from 0.11 to 0.4 fm or so) give a far too large condensate 
term and the wrong sign for the x^ term. 

IV. CONCLUSIONS 

Low-lying eigenmodes of the Dirac operator, which we have previously argued in Ref. |^ couple strongly to instantons 
and anti-instantons, dominate the nonperturbative part of short distance correlators for light quark masses in the 
pseudoscalar, scalar, vector and axial vector channels. In the pseudoscalar channel they are responsible for the bulk 
of the observed strong attraction. The zero eigenmodes in the scalar channel contribute a strong repulsion at larger 
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distances. The pseudoscalar plus scalar channel is more attractive on the lattice than in the instanton liquid model. 
The low modes dominate the vector minus axial vector channel, which has no perturbative contribution. The lattice 
calculations quantitatively reproduce the results of instanton liquid models in the sum and difference of vector and 
axial vector channels. The difference between the lattice results and correlation function extracted from tau data 
in the V + A channel can plausibly be attributed to short distance physics. As the quark mass rises, the relative 
importance of the light modes decreases in all observed channels. 

The lattice simulations could be done much better. It seems to me that these simulations are rather more sensitive 
to systematic effects than they are to statistics. Since one is comparing results to free field theory, an overlap action 
with better "kinetic" (dispersion relation, rotational invariance) properties might be a better choice 1 30 1 . Even though 
these correlators are called "short distance," a larger volume is needed to suppress zero mode effects and to deal with 
the weak (power law) falloff of the normalizing free correlators. (A nonchiral action would be cheaper to simulate, 
but would suffer from fatal exceptional configurations at the low values of quark masses needed to study interesting 
physics questions.) Finally, if one wants to look for OPE terms in the correlators it will be necessary to drop the 
lattice spacing and convert to links which smear only a minimal distance-the hypercubic blocking of Hasenfratz and 
Knechtli would be a good choice. It seems to me since most of the physics which is being explored here involves 
chiral symmetry, it would be a mistake to look at these correlators with a nonchiral lattice action. 

With the use of a chiral fermion action, lattice calculations may be poised to say more interesting things about 
spectral functions in QCD. 
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